# 方程组形式：
# { x₁ + 2x₂ + 3x₃ = 5 }
# { x₁ + 6x₂ + 7x₃ = 9 } 
# { x₁ + 10x₂ + 6x₃ = 8 }

import numpy as np

# 系数矩阵
A = np.array([[1, 2, 3],
              [1, 6, 7],
              [1, 10, 6]])

# 常数向量
b = np.array([5, 9, 8])

# 求解方程组
x = np.linalg.solve(A, b)

print("方程组的解为：")
print(x)
